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We show that the metastable, symmetry-breaking ground states of quantum many-body Hamiltonians have 
vanishing quantum mutual information between macroscopically separated regions, and are thus the most classi¬ 
cal ones among all possible quantum ground states. This statement is obvious only when the symmetry-breaking 
ground states are simple product states, e.g. at the factorization point. On the other hand, symmetry-breaking 
states are in general entangled along the entire ordered phase, and to show that they actually feature the least 
macroscopic correlations compared to their symmetric superpositions is highly non trivial. We prove this result 
in general, by considering the quantum mutual information based on the 2—Renyi entanglement entropy and 
using a locality result stemming from quasi-adiabatic continuation. Moreover, in the paradigmatic case of the 
exactly solvable one-dimensional quantum XY model, we further verify the general result by considering also 
the quantum mutual information based on the von Neumann entanglement entropy. 
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I. INTRODUCTION 


The emergence of a macroscopic classical behavior from 
a microscopic quantum world can be explained in terms of 
decoherence to the environment that quickly destroys the co¬ 
herent superpositions of macroscopic objects (Schrodinger 
cats) m. The selected pointer states must then be factor¬ 
ized states with respect to a tensor product structure that is 
local in real space mis. Similarly, superselection induced by 
decoherence due to weak interactions with the environment 
plays a key role also in the phenomenon of spontaneous sym¬ 
metry breaking, where different ordered sectors with broken 
symmetry are dynamically disconnected and are thus the only 
states that are metastable [41 whence the notion of sponta¬ 
neous symmetry breaking 0. 

In the paradigmatic case of the quantum Ising model, the 
ground space of the ferromagnetic phase at zero transverse 
field h is spanned by two orthogonal product states \0)® N 
and \1)® N which are in the same class of pointer states of 
the typical decoherence argument, while the symmetric states 
\k± = 1/V2(\0)® N =b |l) 0iV ) realize macroscopic coherent 
superpositions that are not stable under decoherence mu. 
Therefore, at zero transverse field h, the situation is very clear: 
the only stable states are those that maximally break the sym¬ 
metry of the Hamiltonian, and at the same time, those that fea¬ 
ture vanishing macroscopic total correlations, including en¬ 
tanglement, between spatially separated regions. 

As we turn on the external field h , we have a whole range 
of values where, before a critical value h = h c is reached, 
there is a magnetic order associated to spontaneous symmetry 
breaking (6|, and the decoherence argument applies within the 
entire ordered phase. This means that, again, the only stable 
states are those that maximally break the Hamiltonian symme¬ 
try mm. However, now the symmetry-breaking states are en¬ 
tangled, and their mixed-state reductions on arbitrary subsys- 
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FIG. 1: A many-body quantum system is partitioned in three distin¬ 
guishable subsystems, A, B, and the remainder E, so that the total 
Hilbert space acquires the tensor product structure PL = PLa^PLb 0 
PLe- The quantity d(A, B) is the distance between the two regions 
A, B , and l is the distance defining the new effective support after an 
adiabatic deformation of operators with initial support on PL a, PL b- 


terns possess in general nonvanishing entanglement (MU, 
as well as quantum [12-14] and classical correlations (6). In- 
deed, the symmetry-breaking ground states can be, “locally”, 
more entangled than some nearby symmetric states f!5l . On 
the other hand, it is always implicitly assumed that such states 
are not macroscopically correlated, while their symmetric su¬ 
perpositions are, in complete analogy with the case h = 0. 
Although this is a very plausible picture, a rigorous proof 
has never been provided, due to the mathematical difficul¬ 
ties in dealing with measures of entanglement and correla¬ 
tions based on the von Neumann entropy; see, e.g., the dif¬ 
ficulties in proving the boundary (area) law in generic gapped 
systems [16. JJ|, or in proving the stability of topological en¬ 
tanglement entropy in topologically ordered states lHH . The 
symmetry-breaking states obey the boundary law for entan¬ 
glement MM, while the macroscopic correlations featured 
by the superposition of two different symmetry broken sectors 
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are of order one. 

The question is then about which quantity one should look 
at in order to distinguish the presence of macroscopic corre¬ 
lations, among all possible source of entanglement and corre¬ 
lations? Historically, the key concepts that have been consid¬ 
ered are the off diagonal long range order (ODLRO) [23 ] and, 
more recently, the two-site concurrence (entanglement of for¬ 
mation) at large distance EUE]. If there is either ODLRO or 
nonvanishing concurrence between two sites in different clus¬ 
ters A and B , then also the two clusters must be entangled, 
since the total state of the global system is pure. This is an 
important point, because the reduced subsystems being in a 
separable state does not imply that there must be no entangle¬ 
ment between the two clusters in a pure state. Even if all the 
remaining correlations are classical, they are due to the fact 
that the overall state is a pure entangled state. 

On the other hand, the reverse argument needs not ap¬ 
ply: it is possible that macroscopic clusters are entangled 
even if measures of two-point correlations, like the concur¬ 
rence, are vanishing. For example, this can happen in the 
2D toric code l26l or in the one dimensional cluster mod¬ 
els l27l - [29l where all two-site concurrences are identically 
zero and yet the macroscopic block entanglement entropy is 
finite. Also states that possess volume law for the bipartite en¬ 
tanglement typically have no two-site entanglement, because 
of monogamy. Moreover, concurrence and ODLRO have to 
be computed for every specific model and cannot provide uni¬ 
versal classifications. 

In the present work, we will show that macroscopic correla¬ 
tions are generally nonvanishing in generic symmetric states, 
and vanishing in symmetry-breaking states. To this end, we 
consider the total correlations (classical plus quantum) be¬ 
tween two generic subsystems A and B of the total system, 
as measured by the mutual information (30j: 

1{A\B) = S(A) + S{B) - S(AB) , (1) 

where S(X) is the von Neumann entropy of the density ma¬ 
trix pertaining to subsystem X , as shown pictorially in Fig.[l] 
The mutual information is indeed a bona fide measure of total 
correlations (classical plus quantum) ED. If for two arbitrary 
subsystems A, B spatially separated by arbitrarily large dis¬ 
tances d(A, B) the mutual information 1(A\B) is vanishing, 
we are assured that there are no macroscopic correlations and, 
in particular, no macroscopic entanglement and no macro¬ 
scopic quantum correlations (all correlation and entanglement 
measures are non-negative defined). Otherwise, taking into 
account that the total system is in a global pure state, the two 
subsystems must be macroscopically entangled and quantum 
correlated. It is immediate to verify that I(A\B) vanishes on 
symmetry-breaking states that are product states of the form 
|0)®* and , while for symmetric superpositions it 

is always of order one, irrespective of the actual value of the 
distance d(A,B). 

We will show that the above result is in fact valid in gen¬ 
eral. In order to prove such statement, our strategy is the 
following. Starting from a factorization point, i.e. a point 
in which the system admits a fully separable pointer state as 
global ground state (32-35 ], we will consider adiabatic de¬ 


formations of the ground state and we will study the behav¬ 
ior of the macroscopic mutual information. The deformation 
corresponds to the adiabatic continuation of the ground state 
obtained by switching on the transverse magnetic field h. We 
will prove that, in the entire symmetry-breaking phase, the to¬ 
tal macroscopic correlations and, a fortiori , the macroscopic 
entanglement, as measured by the mutual information, vanish 
in the maximally symmetry-breaking ground states at large 
spatial separations between arbitrarily selected macroscopic 
subsystems. On the other hand, we will also prove that, as 
long as the deformation is sufficiently small, the macroscopic 
mutual information remains finite in symmetric states. Later 
on, we will apply the results of the general analysis to spe¬ 
cific spin models, showing that the result holding for slightly 
deformed symmetric states is in fact valid for all symmetric 
states in the entire symmetry-breaking (ordered) phase, until 
a quantum phase transition point is reached. 

Proving the general result analytically would be quite a 
daunting task if one were to consider the von Neumann en¬ 
tropy S. Rather, we will resort to a related quantity, the 
2—Renyi entropy, S 2 , and the corresponding 2—Renyi mutual 
information. For specific, computable, examples, we will then 
show that the general conclusions reached using the 2—Renyi 
entropy hold as well using the von Neumann entropy. 

The Renyi entropies of index a are defined as 

S a (A) = (1 - a)~ l log 2 Tr p a A , (2) 

and they are all bona fide entanglement monotones (36). In 
particular, the 2—Renyi entropy S 2 is an experimentally ac¬ 
cessible quantity, since it is the expectation value of the swap 
operator on two copies of a quantum state El, while the von 
Neumann entropy S is not an experimentally friendly observ¬ 
able, because it requires complete state tomography, which is 
essentially impossible on a many-body state. Concrete pro¬ 
posals for measuring S 2 resort on quantum switches (38lL or 
on multiparticle interferometry [ 39]. 

We thus consider the 2—Renyi mutual information defined 
as 

1 2 {A\B) := S 2 {A) + S 2 {B) - S 2 {AB) . (3) 

Unlike the quantum mutual information defined in Eq. 0 . 
the 2—Renyi mutual information may not be positive defined 
for every state in Hilbert space. In order to ensure positivity, 
one can regularize this definition in terms of 2—Renyi rela¬ 
tive entropies, as shown in Refs. (40l Effll . However, on the 
class of states of interest for the present investigation, non- 
regularized 2—Renyi mutual information is always positive 
definite. Moreover, as already mentioned, in the paradig¬ 
matic case of the exactly solvable one-dimensional quantum 
XY model, we will also make direct comparison with the von 
Neumann-based mutual information, finding complete agree¬ 
ment. Finally, the clustering property and other general prop¬ 
erties of the 2—Renyi mutual information that are at the core 
of our investigation and that we will prove below, can be 
proven without too much effort to hold valid also for the reg¬ 
ularized versions. Therefore, in the following, also in order 
to avoid unnecessary formal mathematical complications, we 
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will consider only the non-regularized version of the 2—Renyi 
mutual information. 

By adopting such strategy, we will obtain the following 
main result: in the ordered phase corresponding to the sponta¬ 
neous breaking of some symmetry of a many-body Hamilto¬ 
nian with nonvanishing local order parameter, the mutual in¬ 
formation between two arbitrarily selected subsystems A and 
B separated by a distance d(A,B ) reaches its maximum in 
the symmetry-preserving ground states and is upper bounded 
by exp(— 0(d(A, B))) in the maximally symmetry-breaking 
ground states, i.e. the ground states that maximize the order 
parameter. 

Therefore, we establish rigorously that - at least for some fi¬ 
nite range of values within the phase - spontaneous symmetry 
breaking corresponds to the suppression of macroscopic co¬ 
herent superpositions, and symmetry-breaking ground states 
are the ones selected in real world by environmental deco¬ 
herence, in complete analogy with pointer states. In the fol¬ 
lowing, we specialize to the class of global Z 2 symmetry. In 
particular, we will consider the specific, but paradigmatic, ex¬ 
ample of the of quantum XY models, and we will show that 
macroscopic entanglement survives until a quantum phase 
transition occurs. However, the central elements of this result 
hold in general and the remaining ones can be easily adapted 
to every instance of spontaneous symmetry breaking for dif¬ 
ferent Hamiltonians and different classes of symmetry groups. 


II. CLUSTERING OF 2-RENYI ENTROPY AND 
LONG-RANGE ORDER 

Macroscopic long-range total correlations are revealed by 
a nonvanishing mutual information X between two arbitrary 
regions A and B separated by arbitrarily large distances. If 
X vanishes when A and B are separated by a distance larger 
than what defines the macroscopic interaction scale, then we 
are assured that there is no macroscopic entanglement. 

Let us first consider the case of fully factorized ground 
states. These states are realized at a precise and unique set 
of values of the Hamiltonian parameters in an ordered phase, 
the so-called factorization point or, in spin systems, the fac¬ 
torizing field, first discovered in Ref. l32l . The general theory 
of ground-state factorization in terms of the response to lo¬ 
cal unitary perturbations was fully developed much later, in 
Refs. 133 -351. Factorized ground states can only occur in an 
ordered phase of strongly interacting many-body systems. In¬ 
deed, they are always maximally symmetry-breaking ground 
states with degeneracy equal to the dimension of the symme¬ 
try group and, being fully product states, they have a trivially 
vanishing mutual information X. Remarkably, at the factor¬ 
ization point, all other ground states can always be expressed 
as coherent linear superpositions of the fully factorized and 
maximally symmetry-breaking ground states l42l . Then, by 
construction, symmetric ground states feature a nonvanishing 
X, no matter how large the distance d(A, B) between the A 
and B regions. Moreover, there are no further types of entan¬ 
glement and quantum correlations involved. Therefore, the 
classicality of symmetry-breaking states is immediately veri¬ 


fied at the factorization point. 

As the Hamiltonian parameters are changed adiabatically 
and move away from the factorization point, ground-state en¬ 
tanglement does not come solely from the macroscopic super¬ 
position of disconnected sectors, but also from the fact that 
fully factorized states are no longer ground states. They can 
be represented as U(s)\0 )® N , U(s)\l)® N , where U(s ) is the 
unitary operator that connects the instantaneous ground states 
of the Hamiltonian H(s). The operator U(s ) is an entan¬ 
gling operator, showing that the symmetry breaking ground 
states are now entangled. This raises the problem whether 
they are still only locally entangled or if they have developed 
some macroscopic entanglement. Similarly, for the symmet¬ 
ric states will be U(s)^±, it is not immediately obvious to 
what extent their entanglement is macroscopic or not for a 
generic value of s. In the following, we show that the entan¬ 
glement and the correlations due to the adiabatic continuation 
U ( s ) are local, i.e. they vanish in the limit of arbitrarily large 
spatial separations d(A,B). 

In order to proceed, we need to discriminate clearly be¬ 
tween the macroscopic and the local contributions to X. We 
first recall that by S 2 {A) we mean the quantity S 2 {A) = 
— log 2 Qa where Qa is the purity defined as Qa = Tr [p\], 
and pA is the reduced density matrix from the ground state 
of the total system to the subsystem contained in the finite re¬ 
gion A. To compute the Renyi entropy of order 2 we will 
use the identity Qa := Tr p\ = Ty(SaP (g)2 ), where p 
represents two copies of the original full state on the dou¬ 
bled Hilbert space XL ® XL'. The operator Sa is the permu¬ 
tation operator (swap operator) of order two with support on 
XLaa' only: S A = S A ® I a where S A = ® a eAS a and S a 
is the permutation operator on the a-th spin of the system, 

i*0. S a |il, .. .i a , • • • 5 in) ® \jl 5 • • -ja 5 • • • 5 jn ) — jH3 *• -ja 3 • • • 3 ^n) ® 
|jl, "da-, •••, jn) • 

Next, we exploit a locality result, stemming essentially 
from the Lieb-Robinson bound l43l . Indeed, for the mutual 
information based on the von Neumann entropy, a recent sem¬ 
inal contribution has shown that X(A\B) is an upper bound 
for the two-point correlation functions and a lower bound to 
exponentially decreasing functions of the ratio between the 
d(A,B) and the correlation length [44]. 

Following Hastings and Wen (45), let us consider a many- 
body Hamiltonian sum of local terms, H(s ) = JA hi(s ) with 
a finite gap A E above the low energy sector for some fi¬ 
nite interval of values of the Hamiltonian parameters s (then, 
out of this interval a quantum phase transition may occur). 
Moreover, the local operators are assumed to be bounded: 
||/ii(s)|| < oo. If the ground state of H(s ) is known for 
a particular, fixed set of values of the Hamiltonian parame¬ 
ters, say So, we may obtain it for any other generic set s by 
the quasi-adiabatic continuation U(s ) induced by a continu¬ 
ous deformation of H(s). A local operator Oa with support 
on A transforms as Oa{s ) = (s)OaU(s). The new oper¬ 
ator Oa(s ) has support on the whole Hilbert space. Never¬ 
theless, the locality result implies that we can arbitrarily ap¬ 
proximate it with an operator Oa' {$) that has support only 
over the Hilbert space associated to a region with diameter 
diam(A') = diam(A) + /, as long as l is larger than the cor- 
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relation length £ induced by the gap A E, and by this making 
an error bounded in this way: ||Oa(s) — Oa'(s)\\ < Ke~ 1 ^. 
The constant K grows like l D where D is the spatial dimen¬ 
sion (e.g. the lattice dimension for localized spins). 

Let now p be the ground state of the system for s = so- The 
purity of the restriction of the evolved state p(s) to a spatial 
region C reads 

Qc(s) = Tr [U(s)® 2 p® 2 (U\s))® 2 S c ] 

= Tr [p® 2 (U\s))® 2 S c U(s)® 2 ] 

~ Tr[p® 2 5 c +i(s)] • (4) 

Here Sc+i denotes the permutation operator with support on 
spins that are at most at distance l from C, and the ~ sign 
means that the error is exponentially small in //£. If the sub¬ 
system C is the union of disjoint and macroscopically sepa¬ 
rated subsets, C = A U B, we have 

Qc = Tr[p® 2 (U\s))® 2 S A S B U{s)® 2 ] 

= Tr[p® 2 S A+l (s)S B+ i(s)] • ( 5 ) 

If the distance separating A and B is much larger than l, i.e. 
d(A, B ) >> Z, we can write 

Qc ^ Tr [p® 2 5a+;(s) ®S b +i{s) <g> I E \ , ( 6 ) 

where E is the complement to A and B together, see Fig. [T] 
Assume first that the initial state at s = 80 is one of the com¬ 
pletely factorized ground states (which, we recall, are also 
maximally symmetry-breaking ground states). In this case, 
p(s o) = Pa 0 Pb ® Pe and we obtain 

Qc ~ Tr [p® 2 <S J 4 +;(s)] Tr [p® 2 5 B+ ;(s)] . (7) 


Therefore, the purity in the region C = AUB is the product of 
the purities in the two separated regions A and B from which 
it follows immediately that Z 2 (A|L>) ~ 0. 

Let us next consider the opposite case in which the ground 
state for s = 5o is a macroscopic coherent superposition 
(Schrodinger cat) of the fully factorized ground states. We 
show that any such superposition leads to a non vanish¬ 
ing Z 2 (A|L>. Since we need to consider two copies of the 
ground state, we will label the factorized states by | a) and 
| b), respectively, on each copy. So, we can write | a) = 
\o)a ® \o)b ® \o)e and similarly for | b). Therefore, | p )® 2 = 
J2 a ,b ot a a b \a) A \a) B \a) E ® \b) A \b) B \b) E , where J2 a \a a \ 2 = 
J2b \ a b\ 2 = 1- P° r the sake of an explicit evaluation, con- 
sider the case of a doubly degenerate ground state manifold 
and subsystems A and B with equal size (number of spins). 
For symmetric superpositions we then obtain: 


J 2 (A|B)~log 2 4-21og 2 ^ (a'b'\S A +i{s)\ab) 

aba'b' 

+ lo g 2 Y. {a'b'\S A+ i(s)\ab} 2 

aba'b' 


= log 2 4 + log 2 


Eqbq'b'( a,&/ |^ + K 5 )l a& ) 2 

(T,aba'b'( a ' b '\ S A+l{s)\ab )) 2 


( 8 ) 


Since for s = 0 the expectation values of ( a'b'\SA{Q)\ab) 
are positive definite, by continuity they are still positive for 
a small enough value of s. So we see that, for small s , 
the second term in Eq. [8j is negative but that the T 2 (A\B) 
stays strictly positive, for arbitrary d(A, B). This is also true 
for any non trivial superposition of the symmetry breaking 
sectors given by the amplitudes a a , so that only the maxi¬ 
mally symmetric breaking ones have exactly zero T 2 (A\B) 
(as d(A, B) —>> 00 ). 

In the following, we determine the exact value of X 2 (A\B) 
for macroscopic coherent superpositions with arbitrarily large 
d(A,B) in the entire ordered phase, beyond the perturbative 
case of small 8 in the case in which A and B are made by a 
single spin. We prove this result explicitly for models with Z 2 
symmetry, but the central elements of the proof are valid for 
arbitrary symmetry groups and arbitrary dimension of A and 
B. In fact, our proof implies that for a maximally symme¬ 
try broken ground state all the correlations function between 
two very far subsystems factorize in the product of the expec¬ 
tation value of the local operators. And this implies that the 
mutual information between A and B must vanish. When we 
turn to consider symmetric ground states, some of these lo¬ 
cal expectation values must vanish because the local operator 
does not commute with at least one of the parity operators that 
define the symmetry group of the Hamiltonian. As result the 
mutual information is expected to be different from zero. A 
comprehensive analysis on the dependence on the size of the 
subsystems and/or the symmetry group of the Hamiltonian is 
in progress il46l . 


III. LONG-RANGE MUTUAL INFORMATION IN 
MODELS WITH Z 2 SYMMETRY 


In the following, we focus on spin-1/2 systems with a 
global Z 2 symmetry, thus described by Hamiltonians that 
commute with the parity operator along a fixed spin direc¬ 
tion, i.e. = (S^crf. In such systems spontaneous symme¬ 
try breaking is associated to the presence of a two-fold de¬ 
generate ground state and an off diagonal long-range order 
along a spin direction g\ that is orthogonal to erf. We show 
that throughout the entire ordered phase the long-range mu¬ 
tual information vanishes identically on states that maximally 
break the symmetry, while it remains strictly positive on any 
macroscopic coherent superposition of the two broken sym¬ 
metry sectors. For the sake of simplicity we focus on the case 
in which subsystems A and B are each one made by a single 
spin but the results can be extended straightforwardly to more 
general choices. 

The two orthogonal symmetric ground states \e) and | o) (re¬ 
spectively, the even-symmetric and the odd-symmetric ground 
states) form a convenient basis that allows to write all other 
ground states \g) in the ordered phase as their linear superpo¬ 
sitions: | g(u,v)) = u\e) + v\o). The reduced density matrix 
pc of the two-spin block C = A U B can be expressed in 
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terms of the 2-point correlation functions as follows (9) : 

pdu,v) = \X G lA ' %B (u,v)a l £<j l g , (9) 

lA^B 


where the expectations G lA, ' lB (tz, v)=(g(u , v)\(j l ^<j l ^\g(u^ v)) 

are on products of Pauli matrices and cr^f. As shown in 
Ref. 0, all correlation functions can be associated to spin op¬ 
erators that either commute or anti-commute with the parity 
operator Therefore, the reduced density matrix pc{u,v) 
can be expressed as the sum of a symmetric part that coin¬ 
cides with the density matrix of the symmetric ground state, 
Pc\u,v), and an antisymmetric part, that is a traceless ma¬ 
trix, p^ (u, v). Taking into account the fact that the two sym¬ 
metric ground states fall in two orthogonal eigenspaces of the 
Parity, it is straightforward to verify that p^\u,v) is indepen¬ 
dent of the superposition amplitudes and hence p^\u,v) = 

(s) 

Pc • 

The reduced density matrix of a symmetric ground state 
thus reads 


P ( C = + m n( a A + tJ B)+ m l tJ A a B +™1 °a°b) . (10) 


In Eq. ( [TO] ) ra M is the expectation value of the local oper¬ 
ator that commutes with the parity while m v is the local or¬ 
der parameter. Exploiting Eq. ( fT 0 | ) one can derive the mutual 

(s) 1 -' 

information X 2 } and evaluate its asymptotic expression for 

d(A 1 B ) —y oo! 


4 S) (°c) = log 2 


1 + 


m 


4 

v 


(!+m 2 ) 2 _ 


( 11 ) 


As mn v ^ 0 throughout the entire ordered phase, the above 
relation shows the presence of macroscopic entanglement and 
correlations in the symmetric, coherent superposition ground 
states throughout the entire phase. Up to this point, this result 
is valid for any model with Z 2 symmetry. The actual val¬ 
ues of m v of course depend on the specific model considered. 
Here, we analytically compute the result for the quantum XY 
model. 

The one-dimensional spin-1/2 quantum XY Hamiltonian 
with ferromagnetic nearest-neighbor interactions in a trans¬ 
verse field with periodic boundary conditions reads: 


N 


*=-£ 

i=l 



x 

i +1 


+ 



y 

i +1 


+ ha? 


( 12 ) 

where erf, p = x, y, z, are the Pauli spin-1/2 operators acting 
on site i, 7 is the anisotropy parameter in the xy plane, h is 
the transverse magnetic field along the z direction, and the 
periodic boundary conditions cr ^ +1 =crf ensure invariance of 
the model Hamiltonian under spatial translations. 

Such model can be solved analytically 0 071 ED and, 
hence, the phase diagram can be determined exactly and in 
great detail. In the thermodynamic limit, for any 7 G (0,1], a 
quantum phase transition occurs at the critical value h c = 1 of 


h 
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FIG. 2: (Color Online) Behavior of the mutual information 
(Z^^oo)) between two spins at infinite distance in the symmet¬ 
ric ground state of the one-dimensional XY model (thermodynamic 
limit) as a function of the anisotropy 7 and transverse field h in the 
ferromagnetic phase 0 < h < h c = 1. Within the ferromagnetic 
phase, it is always m v > 0, and hence ( 00 ) > 0. On the other 
hand, m 0 either at h c = 1, or at 7 = 0. Only at these points 

7 ° (00) = 0. 


the transverse field. For h < h c = 1 the system is ferromagnet- 
ically ordered and is characterized by a twofold ground-state 
degeneracy such that the Z 2 parity symmetry under inversions 
along the spin -z direction is broken by some elements of the 
ground space. Using the analytical solution, in Fig. [2] we have 
plotted the behavior of X 2 ^ (00) in the ferromagnetic phase. 

Given the two symmetric ground states, the so-called even 
\e) and odd 1 0 ) states belonging to the two orthogonal sub¬ 
spaces associated to the two possible distinct eigenvalues of 
the parity operator, any symmetry-breaking linear superposi¬ 
tion of the form 


\g{u,v)) =u\e)+v\o) 


(13) 


is also an admissible ground state, with the complex super¬ 
position amplitudes u and v constrained by the normalization 
condition M 2 +M 2 = i. Taking into account that the even and 
odd ground states are orthogonal, the expectation values of op¬ 
erators that commute with the parity operator are independent 
of the superposition amplitudes u and v. On the other hand, 
spin operators that do not commute with the parity may have 
nonvanishing expectation values on such linear combinations 


and hence break the symmetry of the Hamiltonian Eq. (12). 

In the asymptotic macroscopic regime d(A, B ) —>> 00, 
the general two-spin reduced density matrix for an arbitrary 
ground state reads 


pdu,v) = + 1 ( uv* + vu*) [m u ((j u A + a v B ) 

+ m ll m v (<7 l X<7 B + <j\a B )} . 


(14) 


The corresponding expression for the mutual information 
^ 2 ( 00 ) reads 


J 2 (oo) = log 2 


1 


7777(1 — [uv* + vu*y) 

(1 + m 2 + (uv* + vu*) 2 rnl) 2 


(15) 
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FIG. 3: (Color online) Behavior of the mutual information based 
on the 2-Renyi entropy, X 2 (A\B) (left), and the von Neumann en¬ 
tropy, X(A\B) (right), as functions of the logarithm of the inter-spin 
distance r for different superpositions in the ferromagnetic phase of 
the one-dimensional XY model at (7 — h — 0.5). Black circles: 
u — 1, v — 0 (symmetric state); red squares: u — cos( 0 . 057 t), 
v — sin(0.057r); green diamonds: u — cos(0.l7r), v — sin(0.l7r); 
brown up-triangles: u m. cos( 0 . 157 t), v = sin(0.157r); blue down- 
triangles: u — cos( 0 . 27 t), v — sin(0.27r); violet empty circles: 
u — cos( 0 . 257 t), v = sin(0.257r) (maximally symmetry-breaking 
ground state). The two definitions feature the same qualitative be¬ 
havior; in particular, they both vanish in and only in the maximally 
symmetry-breaking ground states. 


Due to the normalization constraint, \u\ 2 + \v\ 2 = 1, the 
fraction in Eq. ( [15] ) is semi-positive defined and vanishes only 
either at m v — 0, i.e. in the disordered, classical paramagnetic 
phase, or when (uv* + vu*) = 1. Therefore, in the ordered 
phase the only ground states with vanishing long-range mu¬ 
tual information, and hence vanishing macroscopic entangle¬ 
ment and correlations, are the maximally symmetry-breaking 
ground states. At the other end of the spectrum, it is easy to 
verify that the maximum of X 2 , for a fixed value of the param¬ 
eters and m u , is always achieved in the totally symmetric 
(even) and antisymmetric (odd) states, the absolute maximum 
being obtained for m^ = 0 and m v = 1. Finally, since the 
one-dimensional XY model allows for the exact evaluation of 
all entropies in the Renyi hierarchy, in Fig. [3] we compare the 
mutual information based on the 2-Renyi and the one based 
on the von Neumann entropy, finding complete qualitative and 
quantitative agreement. In particular, they both vanish in, and 
only in, the maximally symmetry-breaking ground states. Fi¬ 
nally, using the parametrization u = cos 6, v = sin 0, in Fig. [4] 
we also report and compare for completeness the behavior of 
the 2-Renyi based and the von Neumann based mutual infor¬ 
mation as functions of the superposition parameter 0 for the 
two extreme cases of nearest-neighbor distance r — 1 and 
asymptotic distance r = 00, finding again perfect agreement. 


IV. COMPARISON WITH OTHER INDICATORS OF 
MACROSCOPIC COHERENT SUPERPOSITIONS 

Quantum discord (3lJ 09] is a measure of quantum cor¬ 
relations more general than entanglement that may exist in 
mixed quantum states, including separable ones. It is defined 


FIG. 4: (Color online) Behavior of the mutual information based 
on the 2-Renyi entropy, T 2 (A\B) (left), and the von Neumann en¬ 
tropy, X(A\B) (right), as functions of the ground-state superposition 
parameter 0 , with the parametrization u — cos 6, v = sin 0, for 
two extreme values of the distance r between A and B. Dashed 
black curve: r — 1. olid red curve: r = 00 . The absolute mini¬ 
mum is always realized at 6 — 7t/4, corresponding to the maximally 
symmetry-breaking ground state. At this point, both measures of 
mutual information vanish exactly for r = 00 . 


as the difference between mutual information -which accounts 
for all correlations, both classical and quantum- and the opti¬ 
mal classical correlations between A and B , by maximizing 
over all the measurement on B\ Cab = max^} [S(pa) — 

Sc(pAB\{Bk})]- It is therefore sensible to verify whether 
it may be a good quantifier of macroscopic quantum coher¬ 
ence. The long-range pairwise quantum discord between two 
spins in the ground state of the one-dimensional XY chain 
has been recently investigated in Refs. Banana. It turns out 
that such quantity features a long-range behavior quite analo¬ 
gous to that of the mutual information, with the crucial differ¬ 
ence that it vanishes identically in all possible ground states 
as ra M 0. From a mathematical point of view this can 
be easily explained considering that, in such a case, the two- 
spin reduced density matrix in the symmetric ground states 
at asymptotically large inter-spin distance is indistinguishable 
from the one obtained by the symmetry-breaking Gibbs states 
at zero temperature. 

The localizable entanglement in a many-body pure state is 
defined as the maximal amount of pairwise entanglement be¬ 
tween two spins i,j at arbitrary distance that can be achieved, 
on average, by performing generalized measurements on all 
other spins [50, 51]. This naturally defines an entanglement 
length x that diverges in symmetric states. Numerics show 
that in the maximally symmetry-breaking ground states of the 
XY chain the localizable entanglement behaves like the con¬ 
nected spin-spin correlation functions Q xx that are bound to 
decay exponentially at long distance. Therefore long-range 
pairwise entanglement defined via the localizable entangle¬ 
ment features a behavior quite similar to that of the mutual 
information. However, this is just pairwise entanglement (al¬ 
though at long distance), so it does not quite capture the notion 
of a macroscopic superposition. Moreover, the mutual infor¬ 
mation is more readily generalized to thermal mixed states, 
the Gibbs states at finite temperature, for which the evaluation 
of the mutual information presents no particular difficulty. 
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V. CONCLUSIONS AND OUTLOOK 

We investigated macroscopic entanglement [52] through 
the behavior of the quantum mutual information between two 
macroscopically separated blocks of dynamical variables in 
the ground state of many-body systems featuring spontaneous 
symmetry breaking. This quantity detects macroscopic to¬ 
tal correlations, including entanglement. The main result of 
this paper is that in the entire phase with broken symmetry , 
the symmetry-breaking states have vanishing long-distance 
mutual information, while the latter remains finite for any 
non maximally symmetry-breaking superposition, attaining a 
maximum for the totally symmetric states. This fact is easy 
to prove when considering symmetry-breaking states that are 
completely unentangled (fully factorized), whose symmetric 
superpositions are GHZ states. In order to prove this feature 
in the entire ordered phase, a much more challenging task, 
we followed a strategy based on two ingredients: (i) adopting 
measures of mutual information based on the 2—Renyi and 
on the von Neumann entropies, and (ii) exploiting locality re¬ 
sults about quasi adiabatic continuation of quantum states de¬ 
rived by using the Lieb-Robinson bounds l43ll53l 54]. In this 
way we were able to prove that spontaneous symmetry break¬ 
ing selects the many-body states with vanishing long-distance 
mutual information, and thus macroscopically least entangled, 
and therefore most classical. 

In perspective, we are concerned with the investigation 
of several open problems. In particular, it would be in¬ 
teresting to extend our analysis to the case of subsys¬ 
tems of arbitrarily variable size, in order to observe pos¬ 
sible threshold effects, and to generic classes of symmetry 
groups (46) . Moreover, we are interested in studying the 
case of globally mixed states (24) and, in particular, equi¬ 
librium thermal states of models featuring spontaneous sym¬ 


metry breaking below a critical temperature. At thermal 
equilibrium the system will be described by the Gibbs state 
p eq = Z~ x JA a e~P Ei \Ei)(Ei\ a , where with a we have ex¬ 
plicitly labeled different sectors. Below a critical temperature 
T c , if a labels the sectors with broken symmetry, spontaneous 
symmetry breaking means that in every single realization the 
label a will be fixed by the initial conditions. Therefore, our 
statement is that Gibbs states obtained by fixing a feature 
the least long-range mutual information compared to all other 
non-maximally symmetry-breaking Gibbs states [46]. 

From a different perspective, we remark that many-body lo¬ 
calization has recently become a subject of great interest for 
the condensed matter community. Systems with strong disor¬ 
der featuring many-body localization fail to thermalize and to 
obey the eigenstate thermalization hypothesis (55), because 
their eigenstates are more weakly entangled than in typical 
non-integrable systems. It would be interesting to see how 
the techniques developed in the present work might help in 
providing rigorous results regarding the clustering of mutual 
information in such systems. In fact, our techniques can be 
used to investigate also systems that involve long-range en¬ 
tanglement and correlations, such as topologically ordered 
states (56) Ez) and their resilience in the presence of pertur¬ 
bations or at finite temperature [26. 58, 59 ]. 
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